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Amongst the multitude of state reconstruction techniques, 
the so-called "quantum tomography" seems to be the most 
fruitful. In this letter, I will start by developing the mathe- 
matical apparatus of quantum tomography and, later, I will 
explain how it can be applied to various quantum systems. 
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Quantum tomography, like all state reconstruction 
methods, is concerned with the problem of measuring the 
density matrix p of a physical system. The first obser- 
vations on this topic are dated 1957 |Q, but the greatest 
advancements are a result of this decade's work. Also the 
first experiments were performed during this decade . 
The word tomography stems from medicine, as a conse- 
quence of the resemblance between the main formula of 
the first quantum tomography method, i.e. homodyne 
tomography EH! , and the inverse Radon transform used 
in the CAT Nowadays, however, quantum tomog- 

raphy does not have much in common with its medical 
ancestor. Quantum tomography is a general term, re- 
ferring to any state measurement procedure descending 
from an equation of the same form of equation (^). It is a 
versatile and mighty technique, as it can be applied to a 
great variety of systems and as it includes other methods 
as special cases. 

Homodyne tomography became a particular case of 
quantum tomography, when group representation theory 
was employed || . Even though the latter led to a break- 
through in quantum tomography, it was recognized not 
to be the most general approach. There were, in fact, to- 
mographic formulas (formulas like (|])) that could not be 
ascribed to standard group representation theory. Hence, 
in order to include these cases in a general mathematical 
framework, I will introduce some conditions that com- 
prise the definition of group representation. Let, then, Q 
be a group, TL a separable Hilbert space, B(Tt) the alge- 
bra of linear bounded operators defined in all TL, namely 
Banach algebra [0], and T a linear mapping from Q to 
B{TL) (from now on, the word linear, when referring to 
an operator, will always be understood). If we can find an 
irriducible, unitary ray representation D of Q (a ray rep- 
resentation is such that D(g)D(h) — e 1 ^ sh D(g ■ h), with 
£ G K. [0) and six correspondences a, (3 : Q x Q — > K 
and fijTy — > Q, with i,j ~ 1,2, satisfying the equations 

D(g) T(h) = e ia ^ TifM-h- f 12 (g)) , (1) 
D(g) T\h) = e^ T\f 21 {g)-h- f 22 (g)) , (2) 



for every g,h e Q, then we will say that {T(g)} is a to- 
mographic set. If T is an irriducible, unitary ray rep- 
resentation then {T(g)} is a tomographic set (choose 
D = T, fu(g) = g, ■ ■ ■ )■ The converse is not true in 
general, so the requirement of T being a representation 
is more stringent. When TL is finite-dimensional, the hy- 
pothesis that {T(g)} is a tomographic set is sufficient 
to derive (|J), however the case of dim(7i) = oo needs 
a further condition to make sure that every expression 
converges and can be attributed a precise mathematical 
meaning. More explicitly, T needs to fulfill the following 
inequality |llj : 



s 



\(u\T(g)\v)\ 2 < oo V|u») eTL. 



(3) 



Condition (g) could be quite a nuisance, since it must be 
checked for every couple of vectors \u), \v) G TL. Fortu- 
nately, (||) is equivalent to 



3\u),\v)en : IH T G?)H 2 <oo, 



(4) 



as I will demonstrate now. Let us start by admitting 
the validity of for two vectors \w 2 ) (eventually 
coincident) and consider the set V of vectors \v) e TL of 
the form 



\v) = v g D(g) \wi), with ^ \ v g\ < 00 



(•5) 



i being 1 or 2. V is a linear manifold of TL because 
a\vi) + P\v 2 ) belongs to V if \vi), \v 2 ) e V and a, (3 e C, 



since J2„ \ a v ig + P v 2 g \ < ME, \ V H 



oo. The application of any D(h) to a vector \v) G V, 
yields a vector that still belongs to V, in fact D(h)\v) — 
E g v g ^D(h-g)\w i )=Y Jg Vh-i.gj^- 1 ^ D(g)\ Wi ), 
with E ? K-i.ge^M*- 1 -*) | = J2 g \vg\ < oo. The ir- 
riducibility of the operators D(g) implies that V = TL 
or, in other words, that every vector in TL can be written 
in the form (j|). Hence for every \u), \v) e TL, with the 
help of Cauchy's inequality, we obtain 

]T \{u\T(g)\v)f = 
g 

= E E ul l v gi {w 1 \D\g l )T\g)D{g 2 )\w 2 ) 



91,92 



Em) E imt( 5 )mi 2 <oo. 



i 



namely (@) => @ (@ =S> (|) is obvious). 

To prove equation I will work out another identity 
first: 

Assertion 1 

If A is a trace-class operator on TL and {T(g)} is a to- 
mographic set and satisfies (^) (or (^)) then 



£5> (a\T(g)\ Ui ) <«i|Tt(ff)|6), 



(8) 
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(6) 



wii/i fc = (y|T(g)|V>)| 2 indipendent of the choice of 
the normalized vectors \<p), \ip) G TL. 
Proof: Hypothesis (|^) implies 



]T |<a|T( 5 )|&)(c|Tt( 5 )|d)| < 

£ |( a |T(. 9 )|fe)| 2 ^ |<d|T( 5 )|c> 



< oo . 



(7) 



defining unambiguously (a\T(g)\b) (c\T'(g)\d } for 

all |a), |6), |c), \d) G TL, in accordance with [p"i"| . 
Since the complete space 7Y, as a consequence of Riesz- 
Frechet representation theorem ||, is also weakly com- 
plete, we may infer, by virtue of ( fj), that the sequence 
of the partial sums T(g) \u) (v T'(g) is weakly con- 
vergent, as n — > oo, for all \u), \v) G W. The oper- 
ator I uv = J2 g T(g) \u)(v\ T\g) can then be defined 
as the weak limit of such sequence, as n — ► oo. Us- 
ing equation ([!]) and its adjoint and rearranging the 
sum we immediately get D{h) I uv D< (h) = I uv , which, 
due to Shur's first lemma (for a proof in the infinite- 
dimensional case refer to H), is equivalent to I uv — 
k uv I, where k uv G C depends on \u), \ v) and / is the iden- 
tity in TL. Analogously, equation (^) entails that I uv = 
J2 g T\g) \u)(v\ T(g) is a multiple of the identity. Givena 
normalized \w) G TL, k ww may be easily evaluated, in fact 
k ww = (w\k ww I\w) = (w\I ww \w) = J2 g \( w \ T (g)\ w )\ 2 i 
where the series could be interchanged with the in- 
ner product because of the way /„„ is defined. The 
constant k uv may be expressed in terms of k ww in 
the following manner (let \a), \b) be generic vectors): 
(a\I uv \b) = (v\I ba \u) = (v\u) (w\I ba \w) = (v\u) (a\I ww \b) , 
which (v\u) I . The choice 

\u) = \v) = \tp) normalized and the calculation of the 
mean value of the last equation on the normalized vec- 
tor \<p) produces J2 g \{^\ T {g)\^)\ 2 = k w ,w, proving that 
J2 g (vl^Xff)!^)! 2 is indipendent of the vectors \tp), \ip) 
(for as long as their norm is 1) and will therefore be in- 
dicated simply by k. Schmidt decomposition of a trace- 
class operator A, i.e. A = Yli a i\ u i)( v i\> where 
e are orthonormal sequences and J2i a i < 00 > 

a % > Vi, helps showing that J2 g T(g) AT^(g) is mean- 
ingful. It is indeed sufficient to check the absolute con- 
vergence of the expression 



J2(a\T(g)ATHg)\b) 



for all \a), \b) G TL, to insure the validity of the 
definition of '^2 g T(g)AT'(g) as the weak limit of 

YT g T(g)AT*(g), as n->oo. The inequality 

]T |< a |T( 5 )K}( Uj; |rt( 3 )|6)K 

9 

< \(a\T(g)\ Ui )\ 2 J2 \(b\T(9)\vi)f 

a g 

together with a% > and < oo, guarantees that 

the sum of the absolute values of the terms in (||) is 
^ fc^i a i < 00 ■ Because of the absolute convergence 
we can also rearrange the order of the two sums, obtain- 
ing the assertion's thesis: 

I J2 T{9)AT\g) = WT. ^T(g)\ui) (vAT^g) = 
k — k — 

9 * 9 

= ^ a i (vilut) I = TrA. 

i 

Now, finally, equation (^): 
Assertion 2 

The operator identity 



(9) 



holds, when A,T and k are defined as in assertion 1. 
Proof: Let O be an invertible trace-class operator. Using 
(^) twice, it is straightforward to check that 

]T Tr[AT\g)]OT(g) = £ Tr [ T(g)0] T\g)A . (10) 

9 9 

Expanding the trace on the complete orthonormal se- 
quence {\fi)}, with the help of equation (^) again, we 
may write 

1^2 Tr[T(g)0](< Pi \Tt(g)A\<p j ) = 
k — 

9 

= J2{<Ph\OTr[\<Ph)(<Pi\]A\i Pj ) = (<Pi\OA\ipj) . (11) 



Equations © and © give i J2 g Tr[AT^g)}OT(g) = 
OA, which is equivalent to (^|) because of the invertibility 
of O. 

Note [[12| is devoted to a brief comment on some tech- 
nical features of assertions 1 and 2. 

Before we start to examine physical cases, I will cast 
light on some aspects of (|9|). It is well-known that 
Hilbcrt-Schmidt operators form a Hilbert space, usu- 
ally denoted as o~ c (TL), with inner product (A,B) a = 
Tt[A^,B] pSl, and that the space t c (H) of trace-class 
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operators is contained in a c (TL). If we define P(g) = 
T^(g), equation (|^) can formally be rewritten as 

A = Y J (P(9),A) P{g). (12) 

g 

This equality is of simple interpretation: all vectors of 
t c (TL) can be expanded in terms of a closure relation, 
resorting to elements P(g) that are not necessarily in 
a c (H) but that belong to a larger set, in the same way 
that occurs with the expansion of a vector in terms of 
generalized vectors. A comparison will make the situa- 
tion even clearer. If we identify a c (TL) with L 2 (R) (the 
space of square-integrable functions on R), then t c {TL) 
corresponds to the space <S(R) (test functions on R de- 
creasing rapidly at infinity) and B(TL) to the space of 
tempered distributions 5'flR) dual to <S(R). This anal- 
ogy suggests that formula ([L2h is also valid for A 6 <J C (TL) 
if we define the inner product (P(g),A) as the limit of 
(P(g),A n ) , with the sequence of trace-class operators 
{A n } converging to A. Naturally, not all closure relation 
arise from a group. It is possible to use conditions sim- 
ilar to (|l]) and (0) defined on sets that are not groups 
and still obtain (Q). Or even more generally, there are 
cases of spectral decompositions in a c {TL) that do not 
satisfy anything like ([!]) and (Q) (for example the eigen- 
vectors of a self-adjoint operator from a c (TL) to a c (TL) 
give a closure relation in a c (TC)). However, from an op- 
erative point of view, these more general approaches to 
formula (||) are not very useful. Groups are indeed sim- 
ple objects to be dealing with and they produce quite 
a large number of interesting results. Moreover, closure 
relations that do not exhibit a group structure usually 
derive from a group. To clarify this point, we must keep 
in mind that nothing in (|l2] ) guarantees that {P(g)} is a 
complete orthonormal system. We know it is complete, 
because of (|l^) , but we cannot be sure that it is orthonor- 
mal with respect to the inner product ( , ) Q . And, in fact, 
generally it is not orthonormal, i.e. the set {P(g)} is of- 
ten overcomplete. If an orthonormal system is extracted 
from the original overcomplete set, the group properties 
may disappear and we would be left with a spectral de- 
composition that is not generated by a group, even if it 
originated in a group. We will encounter an example of 
this circumstance afterwards. 

Equation (|J) is a closure relation itself. Choosing A — 
\v)(v\, with \v) normalized and arbitrary, it states that 
{P(g)\v)} is a complete set in TL. In other terms, the 
complete set {P(g)\v)} in TL corresponds to the complete 
set {P(g)} in a c (TL). Note that this connection is not a 
consequence of the specific context (Q being a group and 
T satisfying (@) and (|)) in which we proved formulas (H) 
and (g). 

So far, I assumed that Q was discrete. Nonetheless, (j|) 
and (|9|) apply to other situations. It is useful to recall 
that every unitary irriducible representation of a compact 
group is finite-dimensional (meaning that TL is finite 
dimensional) , whereas every unitary irriducible represen- 



tation of a non-compact group is infinite-dimensional 
(with the exception of the trivial representation). Hence, 
for a finite group or a compact Lie group, the mathemat- 
ical problem simplifies: dim(7i) < oo, t c (TL) = cr c (TL) = 
13(H) and the convergence of k is always granted. In 
particular, for a finite group, by tracing both members 
of @, we immediately recognize that k = dil [f^) , with 
[Q] indicating the order of Q . For a compact Lie group, 
formulas (^|) and ([)]) are obtained by substituting rpr J2 g > 
appearing in (|^) and (^|) for a finite group, with J„ dp(g), 
where dp(g) is Haar's invariant measure for Q . Similarly, 
the formal substitution ^ g — * Jg dp(g) allows to write 
(^) and (^J) for a non-compact Lie group, with a discrete 
group as the starting point. A warning is necessary in this 
case, however. For non-compact groups Jg dp(g) is not 
convergent; moreover, differently from the compact case, 
right and left invariance may correspond to two different 
dMflO Hill- Only when they coincide, i.e. only when Q is 
unimodular, formulas (^ and (^) are applicable. 

Equation @ (and (Bh) deserves its own self-existence 
as a pure mathematical result, but my initial goal was 
different. I was concerned with the physical problem of 
measuring the density matrix, and, thus far, (j^) does 
not give us any clue on how to solve such problem. I 
will now show how (^J) is, in actuality, much nearer to 
the solution than what may appear. If TL is the Hilbert 
space associated with the physical system under consid- 
eration, then the density matrix p is an element of t c (TL) 
and, consequently, can be written in place of A in for- 
mula (H). Moreover, if each T(g) is self-adjoint or is a 
function of a self-adjoint operator, then we can evalu- 
ate the trace in (^) over its complete set of eigenvectors 
{|<7,i}}. This operation yields an expression containing 
quantities of the form (g,t\p\g,t), which can be inter- 
preted as the probability that a measurement of T(g) 
gives the eigenvalue t(g) corresponding to \g, t) (the case 
of t(g) degenerate could be treated analogously). These 
quantities are, in principle, experimentally accessible and 
will be indicated with p(g,t). Formula (|9|) then becomes 
P = J2 g ,tP(9,t)[it(g) T Kg)}, where if the group is not 
discrete, or if T(g) has continuous spectrum, sums must 
be replaced by integrals (obviously there could eventually 
be both sums and integrals). The observation that the 
eigenvectors \g,t) and \g',t) are not necessarily different 
even if g ^ g' suggests to divide Q in classes Gi, requir- 
ing the property that all T(g) corresponding to the same 
class have the same eigenvectors. We would then write 

p - E^(M)E seei []K<?)T%)] ^ Z t , t p(ht)K(z,t), 

where the operator K(i,t) is usually called pattern func- 
tion (again, if i is a continuous index then J^. — > J dp(i)). 
The formula 

p=Y,Piht)K(ht) (13) 

i.t 

is the essence of most state reconstruction methods. It 
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states that measuring the probabilities p(i,t) and calcu- 
lating K(i,t) is all that is needed in order to obtain p. 
The peculiarity of quantum tomography is that K(t,t) 
does not need to be determined by solving an inverse 
problem, it is explicitly given by A YlgeQi ^9)^ id)- 

It took quite long to work out ([l3|), but its generality 
will show its strength now that we turn to examples. 

I will begin with the spin case, as it is the least com- 
plex. The (reduced) spin density matrix of one parti- 
cle with spin 5 (integer or half-integer) is defined in 
a Hilbert space H s , with dim(W s ) = 25+1. The 
compact group SU(2) is particularly suited for the case 
of arbitrary 5, since there exists an irriduciblc, uni- 
tary representation of SU (2) in every finite-dimensional 
space. If SU (2) is parametrized with (■#, ip, -0) be- 
longing to [0, 7r] x [0, 2tt) x [0, 2n] and if n is de- 
fined as (cos(/?sin$,sinc/?sin$,cos - $), then the operators 
R(n, ip) = e - "^'™, where S is the spin operator [jl4|, con- 
stitute an irriducible, unitary representation of SU (2) 
and, consequently, a tomographic set in Tig. Formula 
( |l3| ) then becomes 




where df2s is the area element of the unitary spheri- 
cal surface E, Ks(n,m) is the pattern function given 
by (25 + 1) / o 27r d0^1 e iV'(S-S-m) ) and p (ft,m) is the 

probability that m is the result of a measurement of 5 ■ n. 
The calculation of Ks(n, m), the experimental apparatus 
needed to measure p(n, m) and some numerical simula- 
tions can be found in |||l5|]. Because dim(Hs) < oo, 
it is evident that the tomographic set {R(n, ip)} is over- 
complete. It is then possible to choose a finite number 
of operators R(n, ip) and still obtain a closure relation 
in a c (Tis). As previously mentioned, the operation of 
extraction of a smaller complete set from the entire set 
{R(ft, ip)} usually produces a class of operators not cor- 
responding to a group. However, this need not always 
be the case. The dihedral and tetrahedral subgroups of 
51/(2) [|l0| can be represented with a finite number of 
operators R(n,ip), producing tomographic formulas for 
the cases of 5 = \ and 5=1 (for explicit formulas 
and some further considerations refer to ||[l5|). Unfor- 
tunately, there is only a finite number of finite subgroups 
of SU(2); furthermore, the tomographic set associated 
with the tetrahedral group is still overcomplete. There- 
fore, waiving the group structure is a necessity for the 
obtainment of a complete orthonormal system in <7 c (TCs) 
for a generic S. This does not mean that we have to com- 
pletely give up the use of the operators R(n,ip). Know- 
ing in advance that dim(er c (7i)) = (25 + l) 2 , we can just 
choose (25+ l) 2 linearly independent operators R(n, ip), 
that do not need to form or correspond to a group, 
and apply Gram-Schmidt orthonormalization procedure: 



R = R = fl("2,fe)-(Bi,Ji(n2^2)). By 

131 - P^fO^Jl P -° 2 - ll-R(-»2^2)-(Bl,fl("2,V'2))o -Blllo' • ■ ■' 

with ||0|| = yj{0,0) , VO e a c (H). By definition, 
{Bi}, i = 1,2, . . . , (25 + l) 2 , is a basis in a c (H). Nev- 
ertheless also {B\} is a basis in a c (Ti), therefore every A 
in <r c (H) can be decomposed as A = ^2i(B\,A) b\ — 
J2i Tr[Bj A] B\ . Because the operators Bi are linear 
combinations of the operators R{fii,ipi), with a little 
algebra we get A = J2 i Tr[R(n,i,ipi) A] TZi, where the 
(25 + l) 2 operators TZi are linear combinations of the 
operators R'(n,i,ipi) as a result of the reorganization of 
the sum on i. One may check that (TZ\, R(fij,ipj)) = 5ij 
(Kronecker's delta), which means that {TZi} is the dual 
basis of {R(n,i,tpi)} in a c (H). Calculating the last trace 
on the eigenstates \ fti, m), associated with the eigenvalue 
e -W>< m £ ^ ne p era tors R(fii, ipi), with A = p, we attain 
a finite version of (|l4|): 

(2S+1) 2 s 

P= Yl P(nu m) K. s (i,m) , (15) 

i— 1 m— — S 

with K s (i,m) = e~^ m TZ, (the suffix 5 in JC s {i,m) 
is a remainder of the dependence of the operators TZi 
on the dimension of Hs)- Equations having the same 
form of ( p^[ ) or (|l5|) are quite common in the problem of 
spin state reconstruction (for example Jl^-[l^]). Inciden- 
tally, we might also observe that the orthonormalization 
of (25+1) 2 linearly independent projectors \n%, 5)(ri,, 5|, 
i = 1,2,..., (25+ l) 2 , instead of the operators ii(n,, ipi), 
leads to the same results obtained by Amiet and Weigert 

The Hilbert space Tic associated with a system of n 
spins is given by the tensorproduct of the n single-spin 
spaces. If we were to write (g) for the elements of cr c (Ti.c), 
we would tempted to choose Q = SU{2) x . . . x SU(2) 
(n times) and D(g) = T(g) — R(fti, ipi), with 

g E G ( x and denote respectively the direct product of 
groups and the tensor product of operators) . This choice 
would actually give a valid closure relation in a c (Tic) 
(equation (H )), but the corresponding probabilities in 
equation (|13[) would require measurements on single com- 
ponents of the system, which are not always possible. 
This difficulty, at least for systems of spins i, can be 
overcome with a different approach, as illustrated in p5| . 

As a model for systems associated with infinite 
dimensional Hilbert spaces, we can take the space Tio 
of one mode of the electromagnetic field. Although the 
problem is mathematically identical for other systems 
(for example, Tio is isomorphic to the space of a spin- 
less, non-relativistic particle in one dimension), quantum 
optics gives the unique possibility of measuring the equiv- 
alent of linear combinations of position and momentum, 
namely the so-called quadratures = ^(aV^ + ae -1 ^), 
with <j> G M. and a and indicating the annihilation 
and creation operators respectively |2~]] ]. This oppor- 
tunity can be exploited by choosing the non-compact 
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Lie group of translations in the complex plane, with el- 
ements a S C, as Q. Since the displacement operators 
D(a) = exp(a — a* a) form an irriducible, unitary ray 
representation of Q, such that J c d 2 a |(0|D(a)|0)| 2 = it 
(da = d(Rea) d(Ima) is the invariant measure of the 
unimodular group of translations in the plane [|lo| and |0) 
is the vacuum state), we can set D(a) — T(a) — D(a). 
For the purpose of writing (|l3|), we should, however, ex- 
press the tomographic set in terms of quadratures. This 
can be achieved by parameterizing Q with k £ K and 
4> £ [0, 7r), related to a by the equation a = ^fce 1 ^, since 
T(4>,k)=T(a{4>,k)) 
then reads 



Equation (|T^) for this case 



that the mappings D and T are defined on two distinct 
groups, changing ([!]) and (||) appropriately. Apparently, 
this generalization would produce other interesting tomo- 
graphic formulas. For example, the state p of one optical 
mode could be determined by measuring only the pres- 
ence (or the absence) of photons in the state D(a)p D*(a), 
for every a £ C. 

My gratitude goes to Prof. G. M. D'Ariano and L. 
Maccone of the Quantum Optics Group of Pavia. This 
work, without their collaboration, would not have ex- 
isted. 
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dxp(4>, x) K(0, x) 



(16) 



where p(4>, x) is the probability that measuring we get 
x and K(4>, x) = ± / +~ dk JM e ikf - x - x ^ (note that k = tt). 
Equation ([!(]) is the fundamental formula of homodyne 
tomography, which I will not be discussing here, because 
the literature on it is already abundant ( and refer- 
ences therein). 

Homodyne tomography is not the only technique that 
allows to reconstruct the state p of one mode of the 
electromagnetic field. K. Banaszek and K. Wodkiewicz 
showed how p could be determined by measuring, for 
every a G C and n € N, the probability p(a, n) that 
n is the number of photons in the state D(a)pD^ (a) 
|p2| . The same result can be recovered with equation 
(O), if we maintain the choices Q = {Translations in 
the complex plane} and D(a) = D(a), but we select 
the tomographic set {D^(a) e lyaia D(a)}, where y can be 
any real number that is not a multiple of 2tt. It is not 
hard to check that (fi]) and (|^) are satisfied [^3| and that 
/ c d 2 a|(0|Dt(a)e i ^D(a)|0)| 2 = 2{1 _l 03y y Then equa- 
tion (|l^) is simply 



d 2 a 



p(a, n) K y (a, n) 



(17) 



with K y (a, n) 



2(1— cos y) 



Dt(a) 



Differently from (|14[), neither ( pL 6| ) nor (17) can be 
simplified by extracting a complete subset from the en- 
tire tomographic set, since both -Pi (a) = ^"'D^a) and 
P 2 (a) = (^rr^ T ) _ ^Dt(a)e- i 2' ata D(Q!) form orthonor- 



* 2(1— cos y) ' 

mal systems: (Pi (a), Pi (a') ) Q = (Pi(d), Pi(u')) = 
<5(Re a — Re a') i5(Ima- Ima'), with S indicating Dirac's 
delta. 

Some examples have shown how different state recon- 
struction problems can be treated as particular cases of 
a general method, which can be summarized in (|^) and 
in the consequent (|k|). Although the theory developed 
in this letter is quiete comprehensive, there is still space 
for further generalizations. One possibility is to assume 
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